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حکیده 


روش‌های اندکی به امکان‌سنجی استفاده از برونداد مدل‌های پیش‌بینی عددی در تخمین مقادیر از 
دست‌رفته تصاویر سنجش ازدور پرداخته‌اند. بدین منظور در تحقیق حاضر علاوه بر ارزیابی الگوریتم ٩۳۸‏ 
در بازسازی تصاویر امکان استفاده از برونداد مدل پیش‌بینی عددی 1۷15 در تخمین مقادیر مفقود تصاویر 
سنجش‌ازدور بررسی شد. این مطالعه با استفاده از سری زمانی تولیدات ۲ص مودیس در سال‌های ۲۰۰۰ تا 
۰ میلادی و برای منطقه شمال شرق ايران انجام شده است. نتایج شبیه‌سازی‌ها بر اساس شاخص‌های 
اعتبارسنجی 1[۷]58 ۸ و 82 با یکدیگر مقایسه شدند. ارزیابی‌های کمی نشان دادند که روش 5۳۸ با 
مقدار میانگین خطای ۱/۶۸ درجه سلسیوس, 2۱/۹۵ ۷۲5۳( و 1-۰/۷۹ دقت مناسب و عملکرد خوبی 
در تخمین مقادیر مفقود دارد. اعتبارسنجی و مقایسه الگوریتم‌ها در حالت پایه (آزمون ۱) و حالت استفاده 
از برونداد مدل 1۷115 (آزمون ۲) نشان دادند که در صورت نبود تصاویر کمکی مناسب سنجش ازدور 
می‌توان از حروجی مدل ۷1015 در الگوریتم‌های هیبرید و بازسازی تصاویر استفاده نمود. ارزیابی بصری 
تصاویر بازسازی شده نشان داد که اجرای الگوریتم ٩۳۸‏ برای هر دو آزمون, در بافت تصاویر مورد 
مطالعه الگوی مکانی مصنوعی ایجاد نکرد و روند تغییرات مکانی 181 حفظ شد. 

کلیدواژه‌ها: ابرناکی. الگوریتم. داده مفقود. 1.51 ۱/۳15 
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۱۳ جغرافیا و مخاطرات محبطی شماره سی‌ام 


۱. مقدمه 


پایش» تخمین و پیش‌بینی دقیق توزیع مکانی و زمانی پارامترهای هواشناسی و اقلیمی کاربردهای وسیعی در 
مطالعات محیطی دارند. از این میان دمای سطح زمین ("1,51) پارامتر مهمی در کنترل گرمای سطح و تبادل آب با 
اتمسفر است (2013 .21 6 ال). پارامتر مذکور به‌صورت مستقیم و غیرمستقيم در بسیاری از مخاطرات مانند تغییر 
اقلیم» خشکسالی؛ سیل و غیره تأثیر می‌گذارد. در حال حاضر تصاویر سنجشازدور یکی از مهم‌ترین منابع اطلاعاتی 
برای برآورد دمای سطح زمین به شمار می‌روند (2012 ,.21 61 [۳16020). تکنیک سنجش زدور با پایش و سنجش 
پارامترهای محیطی و هواشناسی می‌تواند خلأ مکانی و زمانی موجود در اطلاعات حاصل از روش‌های سنتی 
جمع‌آوری داده را پر کند اما وجود آلاینده‌ها در هوا؛ ابرناکی و از کار افتادن ابزار سنجنده‌ها؛ باعث از دست رفتن 
اطلاعات عظیمی شده و به‌اصطلاح در تصاویر ماهواره‌ای حفره ایجاد می‌شود. تاکنون روش‌های متعددی برای تخمین 
مقادیر از دست‌رفته پيشنهاد شده‌اند. با مرور مطالعات مختلف مشخص شد به‌طورکلی رویکردهای تخمین مقادیر از 
دست‌رفته در تصاویر سنجش‌ازدور به سه دسته مکانی زمانی و يا مکانی- زمانی تقسیم می‌شوند. 

روش‌های زمانی تخمین داده‌های مفقود پالایه‌هایی هستند که بر مبنای محاسبات ریاضی به دست می‌آیند. از 
معروف‌ترین آنها می‌توان به پالایه (1964 ,012۷) عگ 217ز۹۵۷) اشاره کرد. ایشان روش ساده شده برازش حداقل 
مربعات را جهت هموارسازی پیچیدگی‌ها و محاسبه مشتقات مجموعه متوالی از مقادیر پیشنهاد کردند. ,له 61 
(2008) 020 داده‌های ۲1/۸ سنجنده مودیس را با روش پالایه زمانی بازسازی کردند. (2013) ,.21 61 10292007 
به مقایسه روش‌های زمانی پر کردن حفره‌ها در تصاویر سنجش زدور پرداخته‌اند. ایشان در تحقیفات خود به این 
نتیجه رسیدند که کارایی هر تکنیک به کاربرد و طریقه اجرای آنها بستگی دارد. علاوه بر پالایه‌های ساده روش‌های 
پیچبده‌تری مانند موج فوریه (2012 .21 61 30016) و تحلیل هارمونیک (2015 .21 61 210۷) نیز ارائه و روی 
تولیدات مختلف سنجش زدور اجرا شدند. 

دسته دوم یعنی رویکردهای مکانی برای بازسازی پیکسل‌های از دست رفته در تصویر از اطلاعات سایر 
پیکسل‌های همان تصویر استفاده می‌کنند. درون‌یابی ساده مانند نزدیک‌ترین همسایگی؛ درون‌یابی دو خطی 
(121110627» روش اسپیلاین و وزن‌دهی عکس فاصله (12۷۷) از جمله روش‌های مکانی هستند. در اغلب این 
الگوریتم‌ها از ميانگین وزن‌دهی شده استفاده می‌شود (2013 ,.21 61 20۳2). در این میان روش‌های زمین آمار مانند 
کریجینگ نیز توسط (1999) ,۸0۵1 در درون‌یایی داده‌های مفقود مورد استفاده قرار گرفته است. برخی از محققن 
نیز از درون‌یایی طیفی با روش قطعه‌بندی برای تخمین مقادیر از دست‌رفته استفاده کردند (2007 ,.21 61 [[۷12۷6). 


16۳0۵۲۵ 9۲206 0صه1 1 
۸۲۵2۵ 121 2 


سال هشتم ارزیابی الگوریتم و امکان‌سنجی استفاده از ... ۱۳۵ 


در مقایسه با روش‌های درون‌یابی ساده رویکردهایی که از داده‌های کمکی استفاده می‌کنند. بیشتر مورد توجه 
محققین قرار گرفته‌اند. از جمله روش‌های مکانی-زمانی روش تطبیق هیستوگرام کلی (0۳1۷1۱)) و تطبیق هیستوگرام 
خطی محلی ( 11۳1۷ تیم مشترک سازمان ملی هوا فضا/سازمان زمین‌شناسی ایالات متحده (1۹05/۵۸5۸) را 
می‌توان نام برد. این روش‌ها برای پر کردن راه راه شدگی تصاویر لندست پيشنهاد شدند ( تطع۹۵06 > تتعطمه۱۷00 
7 ,۵۷261). به دلیل عملکرد ضعیف روش‌های مذکور در مناطق ناهمگن (20.,)2011 60 60) روش درونیاب 
پیکسل‌های مشابه همسایگی (5۳1) را پيشنهاد دادند. (2012) ,.21 64 20 به کمک تشوری زمین آمار روش بهبود 
یافته ۷5۳ را با عنوان درونیاب پیکسل مشابه همسایگی زمین آماری ]05 ارائه دادند. 

(2014) ,.21 6 بولمصمصعط1۷0 از فرض‌های مشابه روش ]5۳ برای بازسازی تصاویر لندست در شمال ایران 
استفاده کردند. (2016) ,.21 61 067067) الگوریتم مکانی- زمانی مبتنی بر زیرمجموعه‌گیری را برای تخمین مقادیر 
مفقود داده‌های ۷ مودیس پیشنهاد کردند. هر یک از رویکردهای ذکر شده مزایا و معایب خاص خود را دارند. 
عواملی مانند ساختار الگوریتم سناریوهای بکار رفته. نوع متغیر و منطقه مورد مطالعه روی نتایج تاثیرگذار هستند. از 
اینرو بررسی عملکرد الگوریتم‌ها روی متغیرها و مناطق مختلف ضروری به نظر می‌رسد. 

علاوه بر موارد اشاره شده یافتن تصوير مناسب کمکی در روزهای ابری و يا فراهم ساختن سری زمانی طولانی 
به عنوان داده‌های کمکی از جمله مهم‌ترین عوامل تاثیرگذار بر نتایج روش‌های مبتنی بر سری زمانی می‌باشند 
(2013 ,له اه 1620026210 بررسی منابع مختلف نشان داد که در رویکردهای ارائه شده برای حل این مشکل. 
اغلب از تصاویر سری زمانی سنجنده‌های دیگر استفاده می‌شود. اين تصاویر از نظر ساختار طیفی ممکن است با 
سنجنده تصویر هدف متفاوت باشند و بر عملکرد الگوریتم‌ها اثر منفی بگذارند. 

منابع اندکی به استفاده از داده‌هایی مانند خروجی مدل‌های پیش‌بینی عددی در تولید داده‌های پیوسته مکانی - 
زمانی سنجش‌ازدور پرداخته‌اند. مدل‌های مذکور با حل معادلات اولیه حاکم بر جو و زمین به محاسبه پارامترهای 
مختلف محیطی و هواشناسی مانند دما در پیکره شبکه‌بندی شده کره زمین می‌پردازند (2009 ,.21 61 [[06). در این 
راستا محققینی مانند (2004 یله اه عع0ظ) و (2010 ,له 6 1202) از داده‌های مدل پیش‌بینی عددی برای تخمین 
مقادیر از دست‌رفته تصاویر سنجش‌ازدور استفاده کردند. ایشان به این نتیجه رسیدند که در روزهای ابری که یافتن 
تصوير با زمان نزدیک به تصویر هدف مشکل است. استفاده از برونداد مدل‌های پیش‌بینی عددی می‌تواند راهگشا 
باشد. با توجه به مسائل مطرح شده و نیاز روزافزون به داده‌های پیوسته مکانی-زمانی دما و اهمیت آن در تحلیل 
بسیاری از مخاطرات مانند تغیبر اقلیم» حشکسالی و سیل» هدف مطالعه حاضر ارزیابی رویکرد ارائه شده توسط 
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(2016) ,.21 61 167۳067) برای بازسازی تصاویر 1 مودیس و امکان‌سنجی استفاده از برون‌داد مدل ۷۷5 به عنوان 
داده کمکی در تخمین مقادیر از دست‌رفته تصاویر ۲ص مودیس می‌باشد. 

۲ منطقه مورد مطالعه 

طول 0 تا ۱" درجه شرقی و عرض ۲۰ ۳۸ درجه شمالی می‌باشد (شکل ۱ این منطقه با وسعت ۱۳۱۳۰ 
کیلومتر مربع جزو مناطق با توپوگرافی پیچیده و سطوح ناهمگن می‌باشد. بلندترین نقطه منطقه با ارتفاع ۳۳۰۰ متر 
بالاتر از سطح دریا در قله بینالود واقع شده است و خروجی دشت سرخس پست‌ترین نقطه آن با ارتفاع ۰ متر 


۳ 
‌" 
۲۱ 
۹ 
# ۴ 
1 بِ 
قلعت ۳ 
7 ۳ #۳ 
ا ۳ 
۱ ۲ 
۱ 1 
نس 1 4 
9 ا 
۱ 
۳۹ ان 
سس 1 و کته 
0 00 ود .9۶ 9 3 
۹٩۹٩+ #4 7‏ ,هط( 
۱920 ۱440 مه ۱0 21 ۱ 


شکل ۱- موقعیت منطقه مورد مطالعه 


۳ مواد و روش‌ها 
داده‌های مودیس (0۷۲۵۲(۲5۱) 
در تحقیق حاضر از محصولات ۲ مودیس استفاده شد. مودیس با عرض نوار اسکن ۲۳۳۰ کیلومتر و با شرایط 


قرارگیری روی دو ماهوار» تقریباً پوشش کاملی از زمین را در یک روز فراهم می‌آورد (2001 .۵ 64 00حعمصملق). 
این سنجنده دارای ۳۱ باند با قدرت تفکیک طیفی بالا و تفکیک مکانی ۲۵۰ متر (۲ باند)» ۵۰۰ متر (۵ باند) و ۱۰۰۰ 


متر (۲۹ باند) است. ۱۵-۰ درصد از پیکسل‌های تصاویر اصلی مودیس به دلیل شرایط ابری و با ساير دلایل فاقد داده 
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هستند (2014 ,.21 61 02اتع). اطلاعات دمای سطح زمین (181 و قابلیت انتشان یکی از تولیدات سطح ۳ 
مودیس (2۷10111۸2) می‌باشند. تفکیک مکانی داده‌ها یک کیلومتر مربع و بازه زمانی آنها ۸ روزه است. به منظور 
مطالعه روی تصاویر ذکر شده. سری زمانی دمای سطح زمین از سال ۰ ۲۰۱۱۲ از سایت 
( 0962-20۷ 1وم.110://۳00015 .) دانلود شدند. در این تحقیق. مقادیر دما در تصاویر 1ص با استفاده از روش موجود 


در مقاله (2004 .21 66 ۷۷2۲ به شرح زیر از درجه کلوین به درجه سلسیوس تبدیل شدند. 
)۱( 15_- (0]561 + (رگصل< 0610۲ ۷۵) < 1,97۰ 


داده‌های پیش‌بینی مدل 1۷۲۷5 

مدل ۷۲۷۲5[ (نسل پنجم مدل میان‌مقایس عاماا ۱۵۲/۳۲۵ ) برای پیش‌بینی متغیرهای جوی در مقیاس 
سینوپتیکی به کار می‌رود. به‌طور حلاصه مدل 1۷1۷15 یک مدل منطقه محدود غیر ایستا و متأثر از عوارض زمین 
است که برای شبیه‌سازی یا پیش‌بینی گردش میان مقیاس جو طراحی شده است (1995 ,.21 61 [06). در این تحقیق 
داده‌های مربوط به پیش‌بینی مدل 1۷1۷5 از سایت سازمان ۲۲0۸۵ (1005://۷۷۷/۸۷۰۵۵06028.820۷) و سازمان 
کله) ۲۱۲ (بلعتهع۳۸/۵5://۷۷۸۷۷۰/ ) دریافت شدند. اطلاعات دریافت شده تصاویری با تفکیک مکانی ۰/۵«۰/۵ 
دربجهو دوزه زماتی ۲۰۱۰۰-۱۲۸۰۸ بودند که‌بر اسان ظول,و عرضی منطقه مورد مطالفه تهیه شدند به دلیل وخود 
خطاهای قاعده‌مند در برونداد مدل‌های پیش‌بینی عددی» قبل از استفاده از داده‌ها به عنوان ورودی الگوریتم‌هاء بر 
روی آنها پس‌پردازش آماری انجام شد. 

رویکرد بازسازی مقادیر از دست رفته 

بررسی منابع مختلف در زمینه پرکردن حفره‌های تصاویر سنجش‌ازدور نشان داد که هر یک از الگوریتم‌های ارائه 
شده مزایا و معایب خود را دارند. به عنوان مثال روش‌های مکانی مانند درون‌یایی ساده و يا قطعه‌بندی چندمقیاسی در 
عین سادگی اجراء دقت کمی دارند و در مقیاس یک پیکسل قادر به پر کردن حفره نیستند. در عین حال رویکردهای 
پیچیده‌تر مکانی-زمانی مانند 111۷1 011۷1 و روش‌های زمین‌آماری نسبت به روش‌های ساده دقت بالاتری دارند. 
این روش‌ها در برخی موارد قادر به پر کردن حفره بزرگ نبوده و در عارضه‌های خحطی مانند رودخانه دقت بالایی 
ندارند. بررسی‌های بیشتر نشان دادند که روش‌های پیچیده نسبت به تکنیک‌های ساده لزوماً برتری ندارند ( 6 1100 
2 ,۸110500) و روشی برتر است که در عین حفظ دقت مورد انتظان ساده بوده و به زمان و پیچیددگی 
محاسباتی کمتری نیاز داشته باشد. بر این اساس در تحقیق حاضر رویکرد ارائه شده توسط (2016) ,.21 )6 061067 به 
ور ناسا اسان دروکا سا هکس اف ی هواس تست آوووز 


تلد متتمطامومحصنم ق4صه منصوعع0 لقصمتاع 1 
6۵0 ۵۲1۵ظ۵ومصصاخ 1۵۲ مملاه۵۲۵۵1) زاتوتهبتصنا 2 
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استفاده کرده‌اند. این روش برای داده‌هایی با ساختار آرایه چهاربعدی مناسب می‌باشد. در این رویکرد پیش‌بینی مقادیر 
از دست‌رفته در دو مرحله اصلی صورت می‌گیرد. (۱) مرحله زیرمجموعه‌گیری (۲) پیش‌بینی مقادیر از دست‌رفته بر 
اساس زیرمجموعه‌ها. در اینجا توضیحات کلی روش مذکور برای تخمین یک پیکسل فاقد داده شرح داده شده است. 
برای راحتی ارجاع به الگوریتم (2016 .21 64 06506۲) در بخش‌های بعدی متن؛ رویکرد مذکور با عنوان الگوریتم 
پیش‌بینی زیرمجموعه‌ای )٩۳۸(‏ ذکر خواهد شد. 

اگر متغیر (۸۸) (۱ ۲ 6 [ ,5 ,2]2,2 < 2را به عنوان مقدار هر پیکسل در مکان و زمان خاصی در نظر 
بگيريم متغیر (ما 6 ۵ ,وا 6 5ررا ع رو[ > :2] < ۸ آرایه چهاربعدی داده‌ها را با تعداد < و < ۷ < ۷ 
,۷ ۶ و۷ ملفه نشان می‌دهد. در الگوریتم برای پیش‌بینی مقدار مفقود [م20]20,:0,50,2 < 20 در 4 از 
همسایگی‌های مناسب در اطراف 20 به‌صورت ۸ 2 (20) < 8 استفاده می‌شود. در این حالت پیش‌بینی مقدار 
0 < 2 با برازش رگرسیون چارک به مقادیر تصویر کامل و با رتبه‌بندی (۲۳) تصاویر کمکی حاصل می‌شود. 
معادله زیر مدلسازی چارک »0 را نشان می‌دهد: 

۳0( (۵:)۵ ۲ (رم)و < ( 0121۳ 

بر اساس معادله (۲) اگر و2 تخمینی از چارک مقدار مفقود 20 باشد. می‌توان مقادیر همبستگی (0)00م و 
(80) :8 را با حل مسئله کمینه‌سازی تخمین زد(2012 ,۷10[۷۲. بدین ترتیب مقدار پیش‌بینی شده 20 از معادله (۳) به 
مس اد 

20 < 60)۵0( + 706)۵0( (۳) 

در معادله (۳) تمام همسایگی‌های قابل قبول یافت شده ("8) ب‌صورت "رتبه‌بندی شده‌اند. و 70رتبه تصاویر 
با مقادیر مفقود شده 20 می‌باشد. همچنین در رویکرد مذکور چارک مقادیر 0 به‌صورت میانگین مقادیر 
[:۲6 : 2] تخمین زده شدند. 

سناریوهای شبیه‌سازی تصادفی 

در بخش قبل موضوع اصلی الگوریتم ٩۳۸‏ برای تخمین مقدار پیکسل خالی یعنی استفاده از چارک و یافتن 
همسایگی‌های مناسب در تصاویر سری زمانی بیان شد. برای بررسی دقیق‌تر رویکرد مذکور روی تصاویر 1.8 
مودیس به اجرا درآمد. برای این کار از نرم‌افزار ۷۸11.۸3 استفاده شد. ورودی‌های الگوریتم یک تصویر اصلی و 
تعدادی تصویر کمکی بودند. در اين تحقیق دو نوع اطلاعات کمکی برای تخمین مقادیر مفقود استفاده شدند. 
اطلاعات کمکی نوع اول مربوط به سری زمانی تصاویر 1.87 مودیس و اطلاعات کمکی نوع دوم مربوط به برونداد 


اتود )1۵010 ]موطاناه 1[ 
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مدل 1۷۲15 بودند. برونداد مدل 1۷115 با هدف امکان سنجی استفاده از داده‌های پیش‌بینی شده در بازسازی تصاویر 
سنجش‌ازدور در نواحی و تاریخ‌های فاقد داده به کار رفت. پس از آماده‌سازی داده‌ها؛ برای بررسی کارایی رویکرد 
۸ چهار تاریخ انتخاب شدند. هر کدام از این تاریخ‌ها نماینده یکی از فصول سال بودند. در تصویر مربوط به هر 
کدام از آنها؛ ۱۰ ناحيةٌ پیکسل خالی به‌صورت مصنوعی ایجاد شدند. این ناحیه‌ها یک ماتریس مربع شکل در نمایش 
ماتریسی تصاویر هستند که مختصات و اندازه آن‌ها در ماتریس به‌صورت تصادفی انتخاب شدند. حفره‌های مصنوعی 
ب‌صورت تصادفی در منطقه مورد مطالعه به وجود آمدند تا بتوانند محدوده وسیع تغییرات پوشش زمین را در 
مدل‌سازی وارد کنند. شکل (۲) مختصات جغرافیایی یکی از نواحی خالی شده را نشان می‌دهد. در نهایت 
پیکسل‌های خالی هر ناحیه با استفاده از رویکرد ٩۳۸‏ در دو حالت (آزمون (۱): با داده‌های کمکی 141مودیس؛ 


آزمون (۲): با اطلاعات برونداد مدل 0۷185 تخمین زده شدند. 


ارزیابی روش‌ها 

در این تحقیق برای ارزیابی عملکرد رویکرد ٩۳۸‏ در دو حالت مختلف از معیارهای ريشه میانگین مربعات 
خطا ( 0۷15۳ میانگین اختلاف ((۸) و ضریب تعیین (**1) استفاده شد. معیار ]180۷ به جهت تغییرات حساس 
نبوده و به اندازه گیری خطا می‌پردازد. هر چه مقادیر ۷15 به‌دست آمده به صفر نزدیک‌تر باشد. دقت روش 
برآوردی بیشتر است. ضریب تعیین هر چه به یک نزدیک‌تر باشد. دقت روش به کار رفته بالاتر است. مقدار منفی 


کم برآورد و مقدار مثبت آن نیز بیش‌برآورد را نشان می‌دهد. 


)( 7 1 
2 - , ی 
(ز2 4 ۲۷5[ 


۳ (0) 
۸1( < ۰ 


9 (7- 2208 - مایا مر 
(2 - 2۳,02 2(2 - و2۳67 


در روابط بالا ,2 مقادیر پیش‌بینی شده 161 در نقطه و :2 مقدار واقعی آن می‌باشند. ۸۷ تعداد کل داده‌هاست. 


1 1۵0۲ 1۷۲62۵ 500۵۲۵ ۴۲ 
2 ۸۷۰۲۵2۵ 6 
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شکل ۲- موقعیت جغرافیایی حفره مصنوعی ایجاد شده در منطقه مورد مطالعه (ابعاد حفره ۱۵۰*۱۵۰ پیکسل؛ 


مختصات مر کز حفره (۱۰۰۰: 632 در تصویر ۰( ۲ ۱۲*۰ مودیس) 


نتایج و بحث 

در این تحقیق کارایی رویکرد ٩۳۸‏ در تخمین مقادیر مفقود و امکان‌سنجی استفاده از برونداد مدل ۷1۷15 در این 
امر بررسی شد. الگوریتم‌ها به‌صورت دو آزمون با سناریوهای تصادفی بر روی سری زمانی تولیدات 1,8 مودیس 
(۶0۰ تصویر ۸ روزه) در بازه ۲۰۰۰ تا ۲۰۱۰ اجرا شدند. با اجرای الگوریتم‌ها برای هر تصویر ورودی» چندین 
تصوير خروجی حاصل شد که به دلیل حجم نتایح به دست آمده تصویر ورودی و خروجی تاریخ ۱۱,۶ 
برای ارائه نتایج در این بخش درج شده است. 

اجرای رویکرد 9۴۸ با شبیه‌سازی سناریوهای تصادفی 

نتایج اجرای الگوریتم ۸ بر روی یکی از تصاویر سری زمانی 1.5 (تاریخ ۲۰۱۱/۰۵/۱۶) در جدول (۱) ارائه 
شده است. در این جدول, اطلاعات مربوط به آماره‌های خطای (مقادیر متوسط خحطا بیشترین خطا و انحراف معیار 
خطا) هر آزمون به تفکیک ثبت شده است. آزمون (۱) مربوط به اجرای سناریوهای شبیه‌سازی با کمک سری زمانی 
تصاویر 187 مودیس است. این سناریوها برای ارزیابی عملکرد روش ٩۳۸‏ اجرا شده‌اند. آزمون (۲) مربوط به 
اجرای سناریوهای شبیه‌سازی با کمک برونداد مدل ۷1۷15 می‌باشد. این سناریوها به منظور امکان‌سنجی استفاده از 
داده‌های 1۷115 در تخمین مقادیر پیکسل‌های فاقد داده اجرا شده‌اند. با توجه به جدول (۱) نتایج شبیه‌سازی‌ها نشان 
می‌دهند. رویکرد 5۳۸ که در این تحقیق به عنوان روش پایه استفاده شده است دارای دقت مناسبی است. بطوریکه 
مقدار میانگین خطای به دست آمده ۱/2۸۷ درجه سلسیوس می‌باشد. این در حالیست که (2014) ,.1ه 61 دننک 
میانگین خطای موجود در بازسازی تصاویر 181 سال ۲۰۱۱ را ۳0/۲ درجه سلسیوس گزارش کردند و به این نتیجه 


وشییااکک که روش یشان دق ماش با ترلیت خادیهای موه ان قاری ایسان از تصاری او آبازساری نمیا 
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روش 15) 5۵6۸ استفاده کرده بودند. با توجه به جدول (۱) مشاهده می‌شود که مقدار آماره‌های خحطای آزمون (۲) 
بیشتر از آزمون (۱) به دست آمده است. بدین معنی که اجرای الگوریتم ٩۳۸‏ با کمک نقشه‌های برونداد مدل 1۷15 
دقت کمتری نسبت به آزمون ۱ دارد. اين امر می‌تواند به دلیل وجود عدم قطعیت مدل 1۷1۷15 در پیش‌بینی مقدار 
دمای سطح زمین باشد؛ اما به‌طورکلی از نظر بازسازی داده‌های سنجش‌ازدور میانگین خطای الگوریتم‌های آزمون 
شده برای تخمین دما در اغلب پیکسل‌ها در محدوده قابل قبولی قرار داشتند. در منابع مختلف مقدارمیانگین خحطای 
قابل قبول برای دما بین ۱ تا ۸ درجه سلسیوس ذکر شده است (2004 .له 6 2010:7۷20 ,۷۷۵۵۵0 عک «0کناع۳6۲). 


جدول ۱- مقادیر آماره‌های خطای حاصل از اجرای روش ٩۳۸‏ بر روی تصویر 1,81 مودیس (۱۶ می ۲۰۱۱): 


آزمون (۱): میتنی بر سری زمانی مودیس؛ آزمون (۲): مبتنی بر برونداد مدل 1۷۲۲۷۲5 


بیشترین خطا انحراف معیار خطا متوسط خطا خطا (درجه سلسیوس) 
۷/۸۱۲ ۱/۳۳۶ ۱/۸۷ آزمون ۱ 
۱/۰۱ ۱/۹۸۵ ۱۳/۹ آزمون ۲ 


نتایج تصاویر بازسازی شده 1,1 


همانطور که در بخش‌های قبلی ذکر شد به منظور بررسی عملکرد و ارزیابی الگوریتم‌ه اه حفره‌های مصنوعی در 
تصاویر سنجش ازدور ایجاد شدند و سپس مقدار حطای هر الگوریتم با مقایسه مقادیر تخمین زده شده و مقادیر 
واقعی محاسبه شد. شکل ۳- الف تصویر اصلی 1.81 مودیس (تاریخ تصویر: ۱۶ می ۲۰۱۱) و نمونه یکی از 
حفره‌های مصنوعی ایجاد شده در تصویر 1/61 را در حالت واقعی نشان می‌دهد. اندازه حفره ایجاد شده ۵۰۵۰ 
پیکسل می‌باشد. حفره نشان داده شده در شکل. بخش‌هایی از شمال شرقی استان خراسان رضوی را در برگرفته 
است. شکل ۳- ب حالت بزرگنمایی از ناحیه حفره مصنوعی و تغیبرات مکانی 181 درون آن را نشان می‌دهد. 
نتیجه اجرای الگوریتم ٩۳۸‏ نیز با سناریوهای شبیه‌سازی آزمون ۱ و آزمون ۲ به ترتیب درشکل‌های ۳- ج و ۳- د 
نمایش داده شده‌اند. با توجه به شکل ۳-ج و ۳- ده مشخص می‌شود که رویکرد ٩۳۸‏ تمام مقادیر مفقود در تصویر 
را در همه سناریوهای اعمال شده تخمین زده است. بدین معنی که پیکسل فاقد داده درناحیه حفره مصنوعی تصویر 
وجود ندارد. این در حالیست که بسیاری از روش‌های تخمین داده مفقود مانند روش (2004) ,21 6 ۳60 قادر به 
پر کردن تمام پیکسل‌های فاقد داده نیستند. آزمون بصری تصوير در شکل ۳ج و مقایسه آن با حالت واقعی یعنی 
شکل ۳-ب. نشان داد که در برآورد مقادیر مفقود الگوی مکانی روند 7حفظ شده و الگوریتم مذکور الگوی 
ساختگی در تصاویر وارد نکرده است. این الگوریتم به خوبی توانسته در اغلب مکان‌ها با حفظ الگوی مکانی در 


لبه‌های حفره و داخل آن, مقادیر مفقود را بازسازی کند. تنها در قسمت چارک سمت راست و بالای ناحیه حفره 
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الگوی تغییرات دما با تصویر اصلی 1/5۲ (شکل ۳ -ب) متفاوت است. با توجه به شکل ۳-د و مقایسه آن با شکل 
۳-ب چنین به نظر می‌رسد که در آزمون ۲ (حالت استفاده از برونداد مدل) نیز الگوی مکانی تغییرات 181 حفظ 
شده است و حتی در این مورد الگوی تغییرات دما در چارک سمت راست و بالای حفره بهتر از شکل ۳- ج تخمین 
زده شده است. بدین معنی که در این مورد استفاده از برونداد مدل 1۷1۷15 الگوی مکانی روند دما را بهتر از سری 


زمانی سنجش‌ازدور منعکس کرده است. 


(ب) 
۱ (د) 


شکل ۳- نقشه پر شده حفره مصنوعی ایجاد شده در تصویر 1.8۲ (تاریخ تصویر: ۱۶ می ۲۱ با روش ٩۳۸‏ 


100_26300041 ۷ ۷۲ 
-__ 2 


(الف) 


(جج) ..اس اک-9 


(الف) تصویر اصلی و ناحیه حفره روی آن؛ (ب): بزرگ‌نمایی شده ناحیه حفره در حالت واقعی؛ (ج): 
مقادیرشبیه‌سازی شده آزمون ۱؛ (د): مقادیرشبیه‌سازی شده آزمون ۲ 


بررسی توزیع مکانی خطای شبیه‌سازی‌ها 
تصاویر محاسبه و سپس الگوی تغیبرات مکانی درون ناحیه حفره مصنوعی رسم گردید. شکل‌های (۶- الف و ۶- 
ب) توزیع مکانی مقدار خطای درون حفره مصنوعی ایجاد شده در تصویر هآ را نشان می‌دهند. شکل ۶- الف 
مربوط به سناریوهای شبیه‌سازی آزمون ۱ می‌باشد. در این شکل مقادیر حطا نقریباً در محدوده ۰ تا ۷ درحه 


سلسپوس قرار گرفته‌اند. جنانچه دیده می‌شود در اغلب پیکسل‌هاء مقدار خطا در محدوده بین ۰ تا ۲ درجه سلسپوس 
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می‌باشد. همچنین مقدار خطای بیش از 1 درجه سلسیوس در تعداد اندکی از پیکسل‌های مرکز حفره مصنوعی 
مشاهده می‌شود. این مورد ممکن است مربوط به ساختار روش و الگوی زیرمجموعه‌گیری در همسایگی‌ها و یابه 
علت تغییرات شدید توپوگرافی در این ناحیه باشد؛ اما رویکرد ٩۳۸‏ در فواصل کم. الگوی تصویر هدف را به دقت 
بازسازی کرده است. تخمین دقیق مقادیر مفقود مخصوصاً در لبه‌های حفره ایجاد شده همراه با حفظ الگوی روند 
تغییر دما در منطقه از امتیازات روش ٩۳۸‏ می‌باشد که در این مورد با روش پیشنهادی (2013 ,2 21 ع2۵0) مطابقفت 
داشت. حال آنکه چنانچه (2001 ,560 6 ۵0)) در گزارش خود اعلام کزده‌اند» روش مورد استفاده ایشان در 
تخمین مقادیر لبه‌ها دقت کمی دارد. 

شکل ۶-ب تغییرات مکانی مقادیر حطای آزمون ۲ را نان می‌دهد. در این مورد نیز مانند شکل ۶-الف 
بیشترین خطای شبیه‌سازی مربوط به تعدادی از پیکسل‌های ناحیه مرکزی حفره است؛ اما مقدار خطا از حالت آزمون 
۱ بیشتر بوده و بیش از ۱۳ درجه سلسیوس می‌باشد. دلیل این امر می‌تواند استفاده از برونداد مدل ۷1۷15 و 
انتشارخطای حاصل از آن باشد؛ اما با توجه با اینکه محل حداکثر خطا در هر دو آزمون مرکز حفره مصنوعی 
می‌باشد. احتمال می‌رود مشکل ساختاری در الگوریتم 5۳۸ وجود داشته باشد که نیاز به بررسی بیشتری وجود دارد. 
با این وجود در بسیاری از پیکسل‌هاء مقدار خطا اغلب بین ۱ تا ٩‏ درجه سلسیوس تغییر می‌کند که این مقدار بر 
اساس اظهارات (2004 ببلة 6 ۷/۵0: 2010 ,۷۷۵۵۵ 6 ۳6۲8۵1050۳) دقت مناسبی در تخمین مقادیر مفقود دمااست؛ 
بنابراین در مواردی که تصویر مناسب کمکی در سری زمانی 1.۲ یافت نشود. می‌توان از برونداد مدل 1۷1۷15 در 
تخمین پیکسل‌های خالی استفاده کرد. نتیجه به دست آمده در این مورد با نتایج تحقیقات (2004) ,.21 ۵4 «ع13068 و 


(2010) بلق 6 1208 همخوانی داشت. 
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شکل ۶- توزیع مکانی مقدار خطای شبیه‌سازی در پیکسل‌های فاقد داده ناحیه حفره مصنوعی تصویر 1.57 


سنجنده مودیس. شکل الف: آزمون ۱؛ شکل ب: آزمون ۲ 


۱:۶ جغرافیا و مخاطرات محبطی شماره سی‌ام 


اعتبارسنجی و مقایسه شبیه‌سازی‌ها 

در این تحقیق اعتبارسنجی الگوریتم 5۳۸ و د آن با استفاده از شاحص‌های اعتبارسنجی 13۷15۲ ۸۱ و 
۳ انجام شد. مقدار هر شاخص در جدول (۲) نشان داده شده است. معیار 18۷15۳ به منظور ارزیابی احتلاف میان 
مقادیر واقعی تصویر قبل از ایجاد حفره مصنوعی و مقادیر تخمین‌زده شده برای همان پیکسل‌ها توسط الگوریتم‌ها 
استفاده شد. هر چه مقدار 19۷1517 بیشتر باشد. خحطای تخمین بیشتر است. طبق جدول (۲. مقدار 18۷15۳7 آزمون 
یک کمتر از آزمون ۲ است. بدین معنی که دقت الگوریتم در آزمون ۱ بیشتر است. معیار دوم میانگین اختلاف (۸) 
بود که برای ارزیابی بیش‌برآورد و کم برآورد الگوریتم‌ها مورد استفاده قرار گرفت. مقدار مثبت ۸3 بیش‌برآورد و 
مقدار منفی آن کم برآورد را نشان می‌دهد. در هر دو مورد الگوریتم ۸ کم برآورد داشت. البته در آزمون ۲ یعنی 
استفاده از برونداد مدل ۷1۷15 این مقدار بیشتر است. همچنین در این تحقیق از مقادیر "1 برای مقایسه همبستگی 
میان داده تخمین‌زده شده توسط الگوریتم و مقادیر واقعی استفاده شد. هر چه "18 به یک نزدیک‌تر باشد الگوریتم 
عملکرد بهتری دارد و دقیق‌تر است. مقادیر 157 جدول (۲) نشان دادند که اجرای الگوریتم با سناریوهای مربوط به 
آزمون ۱ با مقدار نزدیک ۰/۸ دقیق‌تر از آزمون ۲ است. مقادیر ۳7 به دست آمده در این تحقیق مانند تحقیقات ( 6 
3 ,.21 ۵6 2۵08 :2012 ,له 61 20 :2011 ,.21 61) بود. محققین مذکور نیز از رویکردهای مکانی- زمانی استفاده 
کرده بودند؛ بنابراین با توجه به نتایج به دست آمده. می‌توان روش ٩۳۸‏ را برای تخمین مقادیر مفقود و بازسازی 
تصاویر ىا در منطقه مورد مطالعه به کار برد. 


جدول ۲- نتایج اعتبارسنجی سناریوهای شبیه سازی با دو نوع داده کمکی روی تصاویر دمای سطح زمین 


شاخص اعتبارسنج 
متغیر مورد مطالعه سناریوهای شبیه‌سازی 
۹1۲( ۳۹9 12 
آزمون ۱ ۱/۹۵ ۲ ۷۹ 
دمای سطح زمین ِ 
آزمورن ۷ ۳/۲۰ ۳- 1/۵۹ 


۵ جمع بندی 

در این مطالعه. عملکرد روش مکانی- زمانی بازسازی تصاویر سنجش‌ازدور با نام ٩۳۸‏ که توسط .21 6 06067 
(2016) پیشنهاد شده بود. تحت دو سناریو مورد بررسی قرار گرفت. نتایج اجرای رویکرد ٩۳۸‏ در هر دو حالت 
نشان دادند که اولاً رویکرد انتخاب شده با دقت بالایی قادر به تخمین مقادیر از دست رفته است و انیا مقدار حطای 
تفت اسله در محدوده خطای قابل قبول در تخمین داده‌های دما قرار دارد (2010 ,۷۷۵۵۵ > ۳6۲20907 
همچنین اجرای الگوریتم ۸۵ با کمک نقشه‌های برونداد مدل 1۷*15 دقت کمتری نسبت به آزمون ۱ داشت؛ اما در 
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عین حال مقادیر ميانگین حطا در این حالت نیز قابل قبول بود. از نقاط قوت الگوریتم مذکور می‌توان به حفظ الگوی 
مکانی روند 11 در منطقه و لبه‌های حفره اشاره کرد. اینن در حالیست که بسیاری از رویکردهای پیشین چنین 
قابلیتی ندارند. مقایسه مقادیر شاخص‌های اعتبارسنجی تحقیق حاضر با ساير تحقیقات انجام شده حاکی از دقت 
خوب شبیه‌سازی‌ها در تخمین داده‌های مفقود بودند؛ بنابراین الگوریتم ٩۳۸‏ در موارد مشابه می‌تواند برای منطقه 
مورد مطالعه استفاده شود. همچنین بررسی‌ها نشان دادند که در مناطق با تویوگرافی پیچیده و ناهمگن, مقادیر خطای 
بازسازی افزایش یافت که از نقاط ضعف روش در مناطق ناهمگن به حساب می‌آیند. پیشنهاد می‌شود استفاده از 
روش مذکور در مناطق ناهمگن و حفره‌های بزرگ با احتیاط انجام گیرد. 

با توجه به نتایج به دست آمده در صورت نبود اطلاعات کافی برای دما در یک منطقه می‌توان از داده‌های مدل 
5 برای پر کردن خلاً موجود در تصاویر استفاده کرد. چرا که بسیاری از مطالعات نیازمند حفظ پیوستگی مکانی- 
زمانی داده‌ها است. این مورد بخصوص در روزهای ابری سال که درصد پیکسل‌های فاقد داده زباد است. اهمیت 
بسزایی دارد. به‌طورکلی رویکردهای تخمین مقادیر مفقود بسته به سناریوهای آزمون و روش‌های مقایسه عملکرد 
مختلفی دارند. همچنین عواملی مانند نوع متغیر و منطقه مورد مطالعه نیز در دقت روش‌ها بی تأثیر نیستند؛ بنابراین 
پيشنهاد می‌شود در تحقیقات آینده رویکردهای مکانی- زمانی دیگری نیز در بازسازی داده‌ها بررسی و با رویکرد 
۸ مقایسه شوند. در این مورد ممکن است توجه به برهم کنش اطلاعات جغرافیایی» شرایط اتمسفری و 


هواشناسی به عنوان متغیرهای ثانویه دقت تخمین‌ها را افزایش دهد. 
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6 10۲ اه عطلووهو مامصصع مج مملامل0عتن تمطلهه۳ عصتصتطامهمن .2004 ر.وی رصهوفتامصفه۴ 
۶ لته اهممتهصما مها ما ممتاهء‌نامموه :فعاه ممتاهتانكمهتامممه ۵ مملاهانه له 
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1053-64 
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همه 1290-5011200 مطع ۵۶ صمتاه2۱1 ۷ .2002 ,ما2 بلنا 200 .00 رعط22 .۷ معط22 .2 م2 ۷۷ 
عامصمک ,عاهل ععامصم0 )ما۵0 عمط ممتتتا0 مک میه/۱۷۱۵0۵ 1۵۳۲۵ م0 ۲۵۱۲1۵۷۵0 مامنال 2۳0۵ 
163-0۰ :(83)1 ماصمصصم1/ ۵( ۵۶ همتعممه 

۷۲۵۲5 عمط ۵۶ م۷۵۱102)0 مه ممحرووعووه باتلقتام) ,200۳4 و2 بلنا رم رع202108 .لا م2021 و2 بط2 ۷۷ 
,201-4 :(23)1 ,۹605118 ع)محصمک 0 لقتنم( لقممتامصعا عممم‌ما ممهتاه معا لدماماع 

26182, واملم عمتوعله عمتنع۴۵۵0۷ .2013 وم ,220 200 ,۲۱ بطه‌ط9 بت)‎ 1۵۲ 1۵22052) ۱1+ ٩], 
1۳2۸۵22617 عامصف .0مطصامصظ ممتاه12هاناعع ه همه فلو ررقم مملفومروم لهع0مرصصما- نامر تفن‎ 
روم ۱ص ۵۶ موجه‎ 131: 182-4 

6۰ مصصتا ۱۱۲۳۷۲ دآنا۱۵ هماع ۵۶ ممتامتاتعمهمم‌ع۳ ,2015 و۷ ملاصمصم]۱۷ 2040 با مق1ل و.ل بتا200 
۶ عطلودعگ مامصمع ,(19 آطضص) وعلیمو مصصلا ۵ ولو لقمه متممصصقط ‏ ۵ ۳۵۲۲۵۲۴۵۵۵۵۵ 
217-۰ :163 ,۱۱۱۱۷۲۱۲0۴۵60 

2 1 ومهع عصتالز رم جمهمتممه لمعتادلاهاوهعع 6۷ ۸ .2012 ,.ل مصقطان) 200 و.ما بتالنا ویک بتا2 
49-0۰ :124 راصعصصمم ربص ۵ عمتوصهد عامصصم۴ .عمه1۳2 51/01۴ + ۲(۷۲ 


